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ABSTRACT 


This  thesis  investigates  the  computer  aided  measurement  of  fiber 
diameters  by  laser  diffraction.  The  proposed  system  consists  of  a  light 
sensitive  Random  Access  Memory  (RAM)  chip  which  collects  light  intensity 
data  from  the  laser  diffraction  pattern.  Measurements  of  the  spatial 
location  of  the  nodes  of  the  diffraction  pattern  enables  the  calculation  of 
the  fiber  diameter.  These  measurements  may  be  performed  manually  which 
is  tedious  and  requires  subjective  judgement  of  the  nodes.  The  alternative 
method  of  direct  processing  of  the  intensity  pattern  was  investigated. 
Simulation  is  conducted  to  examine  the  feasibility  of  this  method.  Results 
show  such  a  system  to  be  capable  of  providing  one  order  of  magnitude 
greater  accuracy  than  optical  microscopy  measurements  (with  a  shearing 
eyepiece)  and  double  the  accuracy  of  manual  laser  diffraction  methods 
with  the  added  advantage  of  permitting  the  option  of  total  computer 
automation  in  data  interpretation.  _ __ _ 
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I.  INTRODUCTION 


Fiber  reinforced  composites  are  replacing  structural  metal 
components  in  today's  aircraft.  The  high  strength  and  reduced  weight  of 
composites  results  in  less  drag,  increased  payload,  and  longer  fatigue  life. 
Additionally,  the  directional  properties  of  composite  materials  provide 
unique  design  advantages  over  conventional  materials. 

As  with  many  developing  technologies,  the  reasons  for  the  success  of 
fiber  composites  was  not  initially  appreciated.  Tsai  [Ref.  1:p.  2]  states 
that  fortunately  the  modem  composite  was  so  strong  it  was  reliable  and 
competitive  in  spite  of  less  than  optimum  design  practice.  Over  the  last 
twenty  years,  much  progress  has  been  made  in  understanding  the  micro  and 
macro  mechanics  of  composite  materials.  As  this  knowledge  matures  and 
is  incorporated  into  the  design  process,  the  full  potential  of  these 
materials  can  be  realized. 

An  important  contribution  to  the  design  process  is  modeling 
structural  reliability.  Development  of  probabilistic  models  for  a 
composite  must  take  into  account  the  complex  relationships  that  exist 
between  fiber  and  matrix. 
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One  important  parameter  in  the  probabilistic  model  is  fiber  diameter, 
since  variations  in  fiber  diameter  affect  fiber  failure  density.  Therefore, 
more  accurate  fiber  diameter  measurement  results  in  an  enhanced 
reliability  model  and  a  more  quantitative  prediction  of  structural 
reliability. 

The  purpose  of  this  research  is  to  investigate  a  computer  aided 
method  of  fiber  diameter  measurement.  The  existing  procedure  of 
diameter  measurement  by  laser  diffraction  is  accurate  to  within  0.5£ 
[Ref.  2:p.  210],  It  is  desired  to  improve  this  accuracy  by  interpreting  the 
diffraction  pattern  with  a  light  sensitive  RAM  chip. 

The  presentation  begins  with  a  discussion  of  diffraction  pattern 
analysis.  This  is  followed  by  an  introduction  to  the  Micron  Eye's  theory 
and  operation.  Next,  the  simulation  is  discussed  in  depth  followed  by  the 
results  and  conclusions. 
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11.  BACKGROUND 


Fiber  diameter  measurement  by  laser  diffraction  is  conducted  by 
obstructing  a  collimated  laser  beam  with  a  fiber  sample.  A  diffraction 
pattern  results  and  is  characterized  by  alternating  maxima  and  minima 
symmetric  about  a  central  maximum.  (See  Figure  1) 


Figure  1.  Fiber  Diameter  Measurement  by  Laser  Diffraction 


In  previous  work.  Perry,  Ineichen,  and  Eliasson,  [Ref.  3]  and  Bennett, 
[Ref.  4!,  interpretation  of  the  diffraction  pattern  consisted  of  finding  the 
distance  between  interference  nodes  (minima).  This  distance  is  related  to 
the  fiber  diameter. 
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Bennett  used  the  slit  approximation  (Appendix  A)  to  relate  the 
distance  between  nodes  to  the  fiber  diameter.  Perry,  et  al.,  introduced  a 
more  exact  solution  by  Kerker,  [Ref  5:p.  260]  and  compared  it  to  the  slit 
approximation.  Perry,  et  al.,  concluded  that  the  slit  approximation  should 
be  treated  with  caution.  Therefore,  Kerker's  solution  has  been  adopted  for 
this  study. 

In  his  paper,  Bennett  successfully  demonstrated  the  feasibility  of 
using  the  light  sensitive  RAM  chip  for  diameter  measurements  [Ref  4].  An 
IS32  OPTICRAM  MICRON  EYE,  manufactured  by  Micron  Technology,  Inc., 
was  connected  os  a  peripheral  device  to  an  Apple  II*  computer.  By 
positioning  the  Micron  Eye  so  that  two  interference  nodes  fit  on  its 
surface  (see  Figure  2),  an  “exposure"  of  the  diffraction  pattern  could  be 
printed  by  the  computer  (see  Figure  3).  Analysis  of  the  printed  diffraction 
pattern  gave  the  diameter  of  the  fiber  using  the  slit  approximation  of 
Fraunhofer  diffraction  theory  (see  Appendix  A). 


Figure  2.  Micron  Eye  Array  with  Two  Interference  Nodes 


1 


II!.  DIFFRACTION  PATTERN  ANALYSIS 

Measuring  fiber  diameters  by  laser  diffraction,  requires  an 
understanding  of  the  diffraction  pattern.  This  chapter  introduces  some 
features  of  diffraction  patterns  which  will  be  useful  in  later  analysis. 

A.  FUNCTIONAL  DESCRIPTION  OF  THE  DIFFRACTION  PATTERN 
The  diffraction  pattern  of  a  fiber  can  be  described  as  follows: 

l/l#=  (2/K#Lif)  |  b#  ♦  2  3  bB  cos(nO)  p  (1) 

where  0  =  the  scattering  angle 

K0  *  2if  /A  ( A  =  laser  wavelength) 

b,  =  J#(or)  /  H,G)(a) 

and  a=  ifd/A  (d,  =  fiber  diameter) 

JJor)  are  Bessel  functions  of  the  first  kind, 

w 

are  Hankel  functions  of  the  second  kind. 
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A  formal  introduction  of  equation  (1)  and  the  related  Fraunhofer 
diffraction  theory  is  presented  in  Appendix  A. 

Figure  4  depicts  the  three  dimensional  nature  of  the  diffraction 
pattern.  Equation  (1)  describes  the  intensities  along  one  linear  position, 
or  slice,  of  the  three  dimension  pattern. 


I 


* 


Figure  4.  Three  Dimensional  Diffraction  Pattern 


A  plot  of  equation  (1)  is  the  Intensity  Profile  which  shows  the 


Intensity  Ratio  (!/lt)  versus  the  angle  theta,  for  a  given  fiber  diameter 

and  a  given  fiber  to  screen  distance  L.  Figure  5  shows  a  typical  Intensity 
Profile  for  a  fiber  diameter  of  6.5  microns. 


INTENSITY  RATIO  *10 

0.0  5.0  10.0  15.0  20.0  25.0 


INTENSITY  PROFILE 

TYPICAL 
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B.  FEATURES  OF  THE  INTENSITY  PROFILE  CURVES 


The  intensity  profile  has  several  features  worth  noting.  The  first  is 
the  central  (or  zeroth  order)  maximum.  The  central  maximum  is  the 
largest  peek  in  Figure  5  and  it  is  many  times  more  intense  than  the 
subsequent  maxima.  It  should  also  be  obvious  that  the  profile  is 
symmetric  about  the  central  maximum. 

Other  features  of  the  curves  ore  shown  in  Figure  6.  Figure  6  is  an 
enlarged  view  of  one  of  the  higher  order  maxima  and  associated  minima. 
As  discussed  in  Appendix  A,  the  higher  order  maxima  ore  not  centrally 
located  between  the  minima,  rather  they  are  displaced  slightly  towards 
the  central  maximum.  Therefore,  the  higher  order  peaks  are  asymmetric 
about  their  maxima.  This  asymmetry  means  the  maximum  derivative  on 
the  upslcpe  side  of  the  curve  will  be  greater  than  the  maximum  derivative 
on  the  downslope  side  of  the  curve.  This  difference  is  dealt  with  later 
when  tuning  for  the  optimum  exposure  is  discussed. 

The  minima  in  Figure  6  are  also  distinctive  features  of  the  curves.  In 
Appendix  A,  it  is  shown  that  the  location  of  the  minimo  are  explicitly 
related  to  the  diameter  of  the  fiber  using  the  slit  approximation: 

m  =  *  of  interference  node 

=  sin  I  mA/d  3  A  =  laser  wavelength  (2) 

d  =  diameter  of  fiber 


19 


INTENSITY  RATIO  *10' 
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Figure  6.  Features  of  the  Intensity  Profile  Curve 
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Equation  (2)  is  only  an  approximation  for  a  fiber,  but  it  is  useful  for 
preliminary  analysis.  Such  analysis  includes: 

(1)  determining  where  to  place  the  Micron  Eye  for  data  collection  and, 

(2)  providing  and  initial  guess  of  fiber  diameter  from  diffraction 
pattern  data. 

C.  EFFECT  OF  DIAMETER  ON  THE  INTENSITY  PROFILE  CURVES 

Figure  7  shows  the  effect  of  diameter  on  the  Intensity  Profiles.  Note 
that  successive  maxima  and  minima  ore  further  from  the  central  maximum 
for  fibers  of  smaller  diameter.  This  behavior  will  be  important  to  later 
analysis. 

D.  THE  MAXIMUM  DERIVATIVE  AND  THE  MINIMUM  NUMBER  OF  POINTS 

% 

The  goal  of  automated  fiber  diameter  measurement  is  to  find  the 
diameter  rapidly  and  accurately.  Speed  and  accuracy  can  be  exclusive.  One 
can  imagine  that  collecting  all  the  points  on  the  Intensity  Profile  curve 
will  result  in  an  almost  exact  determination  of  fiber  diameter,  at  the 
expense  of  time.  Two  questions  must  be  answered.  First,  what  is  the 
minimum  number  of  points  to  uniquely  describe  an  intensity  profile? 
Second,  where  on  the  curve  should  these  points  be  obtained? 


o  oc  ~  --  oUNtt 


Figure  V.  Effect  of  Diameter  on  intensity  Profiles  (3-D) 


22 


1.  Minimum  Number  of  Points 


Equation  (t)  consists  of  an  infinite  series.  Determining  the 
minimum  number  of  points  to  uniquely  define  an  infinite  series  is  not  a 
straightforward  procedure.  Cursory  inspection  of  Figure  8  leads  to  the 
conclusion  that  one  point  will  not  uniquely  define  the  curve,  since 
intensity  curves  for  many  diameters  pass  through  the  same  point  (e.g.,  at  0 
=  .04  radians).  If  the  absolute  intensity  is  not  known,  2  points  will  not 
provide  a  unique  solution  either.  For  the  case  of  non-absolute  intensity 
measurements  (as  is  the  case  in  this  experimental  set-up)  3  points 
appears  to  uniquely  define  a  curve. 

The  curve  in  Figure  9  illustrates  this  observation.  An  imaginary 
line  is  fixed  through  three  points  (i.e.,  the  experimental  measurements). 
The  absolute  intensities  of  these  three  points  are  not  known,  and  neither 
are  the  absolute  spatial  locations  with  .respect  to  the  central  maximum 
(i.e.,  0).  Only  the  relative  spatial  locations  among  these  points  are 
known.  Iteration  of  the  diameter  is  equivalent  to  moving  this  line  (with 
the  x  marks  fixed  relative  to  each  other)  along  different  locations  on  the 
curve,  trying  to  match  all  three  points.  This  iteration  can  be  repeated  on 
curves  of  different  diameters  (Figures  10  and  11)  and  no  other  match  can 
be  found.  Thus,  this  argument  defines  three  points  as  the  minimum 
required  to  uniquely  determine  the  intensity  profile  curve. 
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INTENSITY  RATIO  *10' 
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INTENSITY  PROFILE 


Figure  8.  Effect  of  Diometer  on  Intensitg  Profiles  (2-D) 
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2.  Maximum  derivative 


Having  heuristicolly  established  that  three  points  ore  required,  the 
next  step  is  to  determine  the  optimum  location  from  which  to  select  the 
points.  For  example,  if  points  are  selected  at  the  mimima  (or  maxima)  a 
wide  variation  in  theta  results  from  a  small  change  in  Intensity  Ratio. 
This  variation  is  defined  by  the  derivative  of  the  Intensity  Ratio  with 
respect  to  Theta,  and  at  the  extrema  this  derivative  is  very  small.  It  can 
be  shown  that  the  points  on  the  curve  where  the  derivative  is  a  maximum 
will  have  the  least  variation  error.  Therefore,  it  is  desirable  to  use  the 
Intensity  Ratios  corresponding  to  the  maximum  derivatives  as  the 
optimum  intensity  ratios  for  the  three  points.  These  Intensity  Ratios  will 
be  referred  to  as  the  Threshold  Intensity  Ratios. 

Figure  12  shows  a  portion  of  an  Intensity  Profile  with  a  derivative 
curve  (absolute  value)  superimposed.  Vertical  lines  drawn  through  the 
derivative  curve  maxima  intersect  the  Intensity  Profile  showing  the 
location  of  the  optimum  (threshold)  intensity  points. 


INTENSITY  PROFILE  (PERFECT  DATA) 


Figure  12.  Intensity  Profile  with  Superimposed  Deri  votive  Curve 
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The  MICRON  EVE  image  sensor  is  an  optically  sensitive  Random  Access 
Memory  (RAM)  chip  capable  of  sensing  an  image  and  translating  it  to 
digital  computer  compatible  signals.  The  Micron  Eye  was  selected  to 
introduce  automation  to  the  existing  techniques  of  fiber  diameter 
measurement  by  laser  diffraction.  Automation  is  expected  to  increase  the 
speed  and  accuracy  of  diameter  measurements.  This  chapter  introduces 
the  theory  and  operation  of  the  Micron  Eye. 

A.  PHYSICAL  LAYOUT  AND  DIMENSIONS 

The  Micron  Eye  (IS32  QpticRAM  m)  has  two  arrays  each  containing  128 
rows  x  256  columns  of  sensors.  This  application  will  use  only  one  of  the 
arrays.  The  size  of  an  array  is  4420  pm  by  877pm.  (See  Figure  13) 


128  rovs 


Figars  13.  Hicrii  Egc  Array 


Each  sensor  is  a  light  sensitive  element  is  colled  a  pixel.  The 
physical  organization  of  the  pixels  is  shown  in  Figure  14.  The  128  x  256 
elements  actually  map  into  a  129  x  514  ‘cell  placement  grid*.  This 
arrangement  leaves  ‘space  pixels’  in  between  each  pixel  in  the  row 
direction.  The  space  pixels  can  be  set  high,  low,  or  to  the  level  which 
agrees  with  the  maiority  of  its  nearest  neighbors. 

B.  THEORY  OF  OPERATION 

The  pixels  are  capacitors  which  discharge  a  preapplied  voltage  at  a 
rate  proportional  to  both  the  intensity  and  duration  of  the  impinging  light. 
The  voltage  in  on  exposed  capacitor  is  read  and  digitally  compared  to  the 
fixed  threshold  voltage.  If  the  voltage  is  below  threshold  the  pixel  is  read 
as  WHITE.  If  the  voltage  is  above  threshold  the  pixel  is  read  as  BLACK. 
The  digital  comparison  concept  will  be  on  important  part  of  the  simulation 
in  the  next  chapter. 

After  a  pixel  is  read,  the  row  containing  that  pixel  is  refreshed. 
Refreshing  sets  all  pixels  which  are  below  threshold  to  0  volts,  and  all 
pixels  which  are  above  threshold  to  +5  volts.  [Ref.  6] 


29 


IS32  OpticRAM™  topological  information 


Figure  14.  Micron  Ege  Physical  Organization 
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Operation  of  the  Micron  Ege  is  simple.  The  current  research 
configuration  uses  an  Apple  Macintosh  (512K)  computer.  Control  software 
is  provided  by  Micron  Technology,  Inc.  The  Micron  Eye  is  connected 
through  either  the  modem  or  printer  port.  (See  Figure  15)  The  computer 
also  acts  as  the  power  source  for  the  Micron  Eye. 


Figure  15.  Macintosh  and  Micron  Eye 


2.  Exposure 

An  exposure  of  a  portion  of  the  diffraction  pattern  can  be  made  by 
varying  exposure  times.  A  sample  exposure  is  shown  in  Figure  16.  It  is 
important  to  visualize  that  this  exposure  represents  a  cross  section  of  the 
intensity  profile  curve. 
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Figure  16.  Exposure  and  Its  Relation  to  Intensity  Profile 


Longer  and  shorter  exposure  times  will  vary  the  size  of  the  cross 
section.  Varying  exposure  time  Is  equivalent  to  moving  up  or  down  the 
intensity  ratio  axis  of  the  intensity  profile  curve. 

C.  INTENSITY  MEASUREMENTS  WITH  THE  MICRON  EYE 

Bennett  [Ref.  4]  used  the  Micron  Eye  to  measure  the  distance  between 
interference  nodes.  This  approach  was  compatible  with  the  slit 
approximation  requiring  only  the  theta  locations  of  the  Interference  nodes 
to  give  the  diameter.  In  Kerker*s  equation  [Ref.  5:p.  260],  one  cannot  solve 
explicitly  for  d.  The  solution  requires  knowledge  of  the  Intensity  Ratios 
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and  their  respective  theta  locations.  For  the  Micron  Eye,  Intensity  Ratios 
correspond  to  the  user  controlled  exposure  time. 

1.  Two  Approaches  to  Intensity  Calibration 

The  Micron  Eye  must  be  calibrated  to  a  reference  intensity  level. 
One  method  to  accomplish  this  would  require  two  exposures.  One  at  lt  * 

A  and  another  at  \2  *  A.  Delta  represents  an  unknown  level  above  the 

zero  (or  absolute)  intensity.  (See  Figure  17)  The  difference  between 
these  two  values  eliminates  delta. 

A  second  method  requires  a  fiber  of  known  diameter.  One  exposure  of 
this  “calibration*  fiber  will  allow  determination  of  the  absolute  intensity. 


INTENSITY  PROFILE 


Figure  17.  Calibrating  the  Micron  Eye 
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V.  SldULATim 


This  chapter  introduces  the  simulation  of  the  Micron  Eye  /  Laser 
Di  f f racti  on  di  ameter  measurement  system.  The  si  mul  ati  on  answers  many 
questions  relating  to  the  real  system: 

(1)  where  should  the  Micron  Eye  be  placed? 

(2)  what  exposure  is  best? 

(3)  what  accuracy  can  be  expected? 

(4)  is  the  computer  code  valid? 


The  simulation  combines  diffraction  pattern  analysis  with  Micron  Eye 
operation.  The  worth  of  the  results  will  depend  on  the  accuracy  of  the 
simulation. 


A.  COMPUTER  PROGRAMS 

Several  computer  programs  have  been  written  to  conduct  this 
simulation.  These  are: 

(1)  DAT  AMAKR 

(2)  EXPOSURE 

(3)  DIAFIND 


Appendix  B  discusses  these  programs  in  detail. 
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Briefly: 

(1)  DATAMAKR  is  used  to  produce  Intensity  Ratio  versus  Theta  data 
for  a  given  fiber  diameter  d,  and  screen  to  fiber  distance  L.  The 
user  controls  the  number  of  points  produce  and  their  spacing. 

(2)  EXPOSURE  reads  the  Intensity  Profile  data  generated  by 
DATAilAKR  and  recommends  the  optimum  exposure  based  on  the 
average  of  the  intensity  ratios  corresponding  to  the  maximum 
derivatives.  Exposure  will  also  introduce  random  error  into  the 
data,  as  desired. 

(3)  DIAFIND  uses  the  data  generated  by  EXPOSURE  to  find  the  diameter 
of  the  fiber.  The  programs  begins  with  a  user  input  guess  of  the 
diameter  and  conducts  an  iterative  search  /  comparison  routine  to 
find  the  actual  diameter  of  the  data. 


B.  OVERVIEW 

The  general  approach  in  examining  the  proposed  system  will  be  to  consider 
a  typical  carbon  fiber.  This  fiber  has  been  assigned  the  arbitrary  diameter 
of  7.254jim.  Two  additional  fibers  ±20*  of  the  7.254jim  fiber  are  also 
considered.  These  fibers  define  the  range  of  diameters  from  5.803jim  to 
8.705|im.  See  Figures  18,  19,  and  20  for  Intensity  Profiles  curves  for 
these  diameters. 
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LEGEND 
5.803  MICRONS 


THETA  (RADIANS) 


Figure  18.  Simulation  intensity  Profile  (5.8G3jim) 

SIMULATION  INTENSITY  PROFI 


TIIETA  (RADIANS) 


Figure  19.  Simulation  Intensity  Profile  (7.254(im) 


Figure  20.  Simulation  Intensity  Profile  (8.705jim) 

C.  POSITIONING  OF  THE  MICRON  EYE 

In  order  to  measure  the  diameter  of  the  typical  fiber,  a  method  of 
postitioning  the  Micron  Eye  must  be  defined.  Further,  from  this  some 
position,  it  is  desired  to  also  measure  the  *20%  fibers.  This  requires  that 
one  position  of  the  Eye  permit  fiber  diameter  measurements  over  the 
specified  range  of  fiber  diameters. 


i 


Such  a  position  is  defined  such  that  three  data  points  will  be  obtained 
for  any  fiber  in  the  diameter  range,  as  in  Figure  21.  First,  a  simpler  case 
will  be  discussed  to  introduce  positioning  for  a  given  diameter  fiber. 


12  3  TMa  (rWians) 


MferM  Eg« 

Figure  21.  Exposure  for  Three  Data  Points 
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1.  Positioning  for  a  Single  Fiber 

Suppose  *n  exposure  of  the  first  through  the  second  interference 
nodes  (2nd  maximum)  was  desired  for  the  ?.254pim  fiber.  Figure  22  shows 
the  relationship  between  the  fiber  and  the  Micron  Eye: 


Figure  22.  Micron  Ege  Placement  in  Diffraction  Pattern 

0^  =  tan-1  l  Xj  /  L  1  (3) 

0H1  =  tan"1  [  xH1  /  L  ]  (4) 


.59 


or  tan  0f+i-tan  0s  =  1/L  (xhi  -  xs)  =  Ax/L  (5) 

take  Ax  =  4.4  mm  (width  of  Micron  Eye  array) 

With  Ax  fixed,  one  can  adjust  L  so  that  any  desired  0H1  -  ©f  will 
fit  on  the  Micron  Eye's  array.  Recall  that  0mill  =  sin  I  mA/  d  1 

equation  (2),  which  will  give  the  theta  locations  of  the  first  and  second 
interference  nodes. 

0^-  =  sin  (mA/d)  where  m  =  node  * 

A  =  laser  wavelength  (632.8nm) 
d  =  fiber  diameter 


for?.254ym:  0f  =  .087 1 2  radian s 

0^1  =.  17459  radians  (for  is!) 

with  0^  and  0i+1  fixed,  L  is  defined  from  equation  (5): 


L  is  the  distance  of  the  Micron  Eye  from  the  fiber  in  the 
longitudinal  direction.  Referring  back  to  equation  (3),  the  lateral 
placement  of  the  Eye  is  given  (see  Figure  23): 

x|nm||  =  L  ton(0)  =  4.27  mm 


Figure  23.  Micron  Eye  Position 


Thus,  equations  (3)  through  (5)  give  the  Micron  Eye  position  for  any 
desired  d,  0HI,  0i#  L  and  m. 


This  position  is  not  yet  optimum  for  single  diameter  case.  This 
analysis  has  captured  only  two  interference  nodes  and  the  intervening 
maximum.  Any  exposure  will  yield  information  somewhere  between  the 
two  extremes  shown  in  Figure  24.  The  Micron  Eye  images  correspond  to 
the  shaded  areas  of  the  intensity  profile  curve  on  the  right  side  of  Figure 
24. 


Stort 


&  (rWiaas) 

Figure  24.  Short  and  Long  Exposures 


One  method  of  finding  the  points  would  be  to  search  the  Micron  Eye 
array  data  for  a  major  axis.  The  endpoints  of  this  major  axis  are  the 
points  required  for  the  analysis.  (See  Figure  25) 


Figure  25.  Endpoints  of  Major  Axis 


Recoil  from  diffraction  analysis  (Chapter  3)  that  three  points  are 
required  to  define  the  intensity  profile  curve.  Therefore,  more  of  the 
diffraction  pattern  mu*t  Pe  seen  by  the  Micron  Eye.  A  reduction  in  L  will 
yield  the  third  point.  (See  Figure  26)  L  should  not  be  reduced  any  more 
than  is  necessary,  since  fewer  pixels  are  being  used  to  describe  the  data, 
degrading  the  resolution  of  the  Micron  Eye. 


Major  Axis 

endpoints 

Figure  26.  Three  Points  on  Major  Axis 
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2.  Positioning  the  Micron  Eue  for  o  Range  of  Diameters 

Consider  the  positioning  of  the  Micron  Eye  so  that  three  data  points 
con  he  collected  for  a  range  of  diameters.  A  "window'  for  the  array  must 
be  defined  which  will  ensure  three  points  of  data  for  any  diameter  in  the 
specified  range. 

Examine  the  Relative  Intensity  Profile  plot  in  Figure  27.  0  begins 
at  .075  radians  which  excludes  the  the  central  maximum.  The  central 
maximum  Is  so  Intense,  the  Micron  Eye's  fastest  exposure  cannnot  prevent 
overexposure.  Therefore,  the  central  maximum  is  not  considered  as  a 
location  for  the  window. 

The  window  should  also  be  located  In  a  region  where  the  maximum 
intensity  corresponding  the  the  largest  diameter  fiber  is  close  to  the 
maximum  intensity  of  the  smallest  diameter  fiber.  This  ensures  the 
Threshold  Intensity  Ratio  will  be  more  representative  for  all  diameter  in 
the  range.  This  condition  occurs  at  the  higher  order  nodes  (2  or  greater), 
as  Illustrated  in  Figure  28.  The  plot  of  Intensity  ratio  derivatives  also 
shows  this  condition  in  Figure  29. 
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RELATIVE  INTENSITY  PROFILES 


THETA  (RADIANS) 


Figure  27.  Relative  Intensity  Profiles 


DI/D  (THETA)  *10 

0.0  2.0  4.0  6.0  8.0 


Figure  28.  Region  of  Similar  Maximum  Intensities 
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Figure  29.  Region  of  Similor  Intensity  Ratio  Derivatives 
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To  obtain  three  points,  at  about  1.5  maxima  are  required.  Since 
the  distance  between  nodes  is  the  greatest  for  the  smallest  diameter,  find 
the  0  locations  for  the  smallest  diameter  which  gives  1.5  maxima, 
for  5.803um: 

em*2  =  sin  (21/5.803*-)  =  .216 

=  Sin  (31/5.803m0  =  -321  rads 

Half  of  this  range  is  about  .050  radians  and  1.5  maxima  can  be 

approximated  by  .165 - >  .321  radians.  Admittedly,  this  is  not  a 

sophisticated  method.  The  method  is  somewhat  liberal  and  it  was  later 
discovered  that  a  range  of  .165  to  .300  radians  was  better  than  the  larger 
range.  This  is  because  the  minima  are  not  of  interest  so  the  right  hand 
minimum  at  .321  radians  was  discarded,  and  the  range  reduced  to  .300 
radians.  Only  the  part  of  the  curve  where  the  slope  is  a  maximum  needs  to 
be  “seen"  by  the  array. 

A  quick  check  can  be  made  to  insure  that  the  curve  representing 
the  maximum  diameter  in  the  range  fits  inside  the  window, 
for  8.705 Jim: 

6min2  =  Sin  (2  1  /  8.705)  =  .  1 50  rads 

emin3  =  s1n  (3  1  /  8-705)  =  -2 1 6  rads 

6mi.  =  s:.n  (4 1  /  8.705)  =  .286  rads 


47 


The  next  step  is  to  determine  the  distance  L  at  which  the  Micron 


Eye  will  cover  the©  range  of  .165  — >  .300  radians, 
defines  the  longitudinal  position  and  equation  (3) 
corresponding  lateral  position: 


L  =  Ax  /  [tan(0Hf)  -  tairidj)]  =  30  mm 


Equation  (5) 
defines  the 


xi  =  L  tan(0{)  s  5  mm 


(inboard  position  of  array) 


VI.  SIMULATION:  THE  PERFECT  AND  IMPERFECT  DATA  SETS 


Now  that  a  window  has  been  defined  (0  range)  for  tho  diameter  range, 
it  is  necessary  to  construct  perfect  data  for  the  simulation  diameters:  dp 
dj,  and  dj.  DATAMAKR  generates  this  data  for  a  given  d  and  L,  over  the 
specified  range  of  0.  Additionally,  the  number  of  points  over  the  specified 
range  must  be  chosen.  The  number  of  points  defines  the  0  interval  and 
corresponds  to  the  physical  spacing  of  the  pixels  in  the  Micron  Eye  array 
which  is  on  the  order  of  lOjim.  To  match  this  spacing,  400  points  were 
generated  over  the  interval  .165  to  .300  radians. 

A.  DIGITIZATION  OF  DATA 

The  EXPOSURE  program  is  used  to  digitize  the  perfect  data.  Digitizing 
the  data  is  the  simulative  analog  of  the  Micron  Eye  threshold  voltage 
comparison.  In  the  simulation,  any  intensity  ratios  above  a  certain  ratio 
(call  it  the  threshold  ratio)  are  assigned  a  value  of  1  and  any  intensity 
ratio  below  the  threshold  ratio  are  assigned  a  value  of  0. 
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1.  Determining  the  Threshold  Intensitu  Rati 

Before  EXPOSURE  can  digitize  the  perfect  data,  the  threshold  ratio 
must  be  determined.  EXPOSURE  calls  the  subroutine  DERIV  which 
calculates  the  derivatives  of  the  intensity  profile  data  at  each  theta 
location. 

DERIV  proceeds  to  search  for  all  the  local  maximum  derivatives. 
The  average  of  the  maximum  derivatives  is  returned  to  the  EXPOSURE 
program  and  is  used  as  the  threshold  intensity. 

After  the  data  is  digitized  with  respect  to  the  threshold  intensity, 
the  digital  data  is  searched  to  find  the  theta  locations  where  the  intensity 
changes  from  0  -->  1  or  1  ~ >  0..  The  values  for  theta  at  these  locations 
are  averaged  and  this  average  is  taken  to  be  the  theta  location  where  the 
threshold  intensity  is  located.  The  averaging  of  theta  is  based  on  the  fact 
the  Intensity  Profile  curve  is  approximately  linear  in  the  region  of  the 
Threshold  Intensity  Ratio.  Figure  30  shows  the  general  location  of  the 
Threshold  Intensity  Ratio. 

In  the  simulation,  EXPOSURE  finds  three  average  theta  locations 
corresponding  to  the  threshold  intensity  ratio.  These  data  are  stored  in  a 
data  file  which  is  subsequently  input  to  the  program  DIAFIND.  DIAFIND 
prompts  the  user  to  input  an  initial  guess  of  the  diameter  and  proceeds  to 
find  the  diameter  associated  with  the  EXPOSURE  data. 
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3.  INTRODUCTION  OF  ERROR  INTO  THE  PERFECT  DATA 

EXPOSURE  prompts  the  user  for  an  error  input.  If  error  is  desired; 
EXPOSURE  calls  the  subroutine  RANDOM.  RANDOM  introduces  error  based  on 
the  maximum  intensity  ratio  in  the  perfect  data: 

I  =  L  +  l  I  *  ERROR  *  RND  ] 

fir  o  mx 

where  ERROR  is  the  error  to  be  introduced  and  RND  is  a  random  number 
between  ±1  generated  by  the  NONIMSL  subroutine  RANDU.  RANDOM  returns 
the  now  "imperfect"  intensity  ratios  to  the  EXPOSURE  program. 

C.  DIGITIZATION  PROBLEMS  WITH  THE  IMPERFECT  DATA 

Digitizing  the  imperfect  data  rcsults  in  local  regions  where  the 
intensity  oscillates  between  0  and  1,  as  shown  in  Table  1.  The 
introduction  of  random  error  changes  the  smooth  curve  to  erratic  points. 
Locdly,  variations  above  and  below  the  threshold  intensity  ratio  cause  a 
series  of  0~>1  and  1— >0  oscillations.  (See  Figure  31)  Ideally,  only  one 
local  theta  location  is  to  be  associated  with  the  threshold  intensity  ratio. 

There  are  two  approaches  to  this  problem.  First,  the  average  of  the 
local  theta  locations  can  be  calculated.  This  method  results  in  three 
theta  locations  for  input  into  the  DIAFIND  program  (the  same  as  perfect 
data). 
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TABLE  1.  DIGITIZATION  ERROR 


THETA  LOCATION  INTENSITY  RATIO 


0. 2321624000E  00 
0. 2324997000E  00 
0. 2328373000E  00 
0. 2331748000E  00 
0. 2335124000E  00 
0. 2338498000E  00 
0. 2341873000E  00 
0. 2345249000E  00 
0. 2348623000E  00 
0. 2351998000E  00 
0. 2355374000E  00 
0. 2358747000E  00 
0. 2362123000E  00 
0. 2365499000E  00 
0. 2368872000E  00 
0. 2372248000E  00 
0. 2375622000E  00 
0. 2378997000E  00 
0. 2382373000E  00 
0. 2385748000E  00 
0. 2389124000E  00 
0. 2392498000E  00 
0. 2395873000E  00 
0. 2399249000E  00 
0. 2402623000E  00 
0. 2405998000E  00 
0. 2409374000E  00 
0. 2412747000E  00 
0. 2416123000E  00 
0. 2419497000E  00 


0. 3544308000E-04  1 

0. 3196243000E-04  1 

0. 3160309000E-04  1 

0. 3246925000E-04  1 

0. 3372986000E-04  1 

0. 3785736000E-04  0 
0. 3090432000E-04  1 

0. 3342219000E-04  1 

0. 3387519000E-04  1 

0. 3205373000E-04  1 

0. 3516222000E-04  1 

0. 3860490000E-04  0 

0. 3421140000E-04  1 

0. 3468163000E-04  1 

0. 3580747000E-04  1 

0. 3904779000E-04  0 

0. 3578593000E-04  1 

0. 3912353000E-04  0 

0. 3599435000E-04  1 

0. 3474336000E-04  1 

0. 3538676000E-04  1 

0. 3470230000E-04  1 

0. 3952176000E-04  0 

0. 3923557000E-04  0 
0. 3944175000E-04  0 

0. 3966210000E-04  0 

0. 4409726000E-04  0 
0. 4269459000E-04  0 

0. 4371968000E-04  0 

0. 3904802000E-04  0 


INTENSITY  RATIO  *10 

27.0  29.0  31.0  33.0  35.0 


ERROR  SIMULATION 
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The  second  method  is  to  accept  the  theta  locations  as  they  are  and 
enter  them  into  the  DIAFIND  program.  The  latter  method  will  pass  a 
variable  number  of  data  points  to  DIAFIND. 

Both  approaches  were  tested  (at  5%  error  for  7.254  pm)  and  the  seme 
diameter  was  recovered  by  DIAFIND  in  each  case.  The  only  difference 
between  the  two  methods  is  that  the  averaging  method  runs  one  second 
faster  (out  of  22  seconds  on  the  IBM  360)  than  the  other  method. 

The  averaging  method  was  adopted  for  this  simulation  for  two 
reasons: 

(1)  The  data  input  to  DIAFIND  will  always  be  a  constant  number  of 
points  so  that  DIAFIND  can  be  easily  adapted  to  run  actual  Micron 
Eye  data  without  simulation  related  logic  buried  in  the  code. 

(2)  The  averaging  method  returns  the  same  diameter  as  the  other 
method. 

D.  CHOOSING  THE  CORRECT  EXPOSURE  (TUNING) 

The  optimum  exposure  for  this  simulation  was  calculated  by  the 
program  EXPOSURE.  It  was  an  average  value  of  the  intensity  ratios 
associated  with  the  maximum  derivatives  locations  for  the  7.254pm  fiber 
data.  This  same  exposure  was  used  throughout  the  simulation  for  all  the 
fibers. 

Using  the  same  exposure  introduces  realism  into  the  simulation  since 
one  cannot  tune  the  exposure  for  every  fiber  that  is  to  be  measured. 
Because  the  Micron  Eye  window  is  located  at  least  two  nodes  out  from  the 


central  maximum  the  exposure  is  closer  to  optimum  for  all  diameters  over 
the  range.  (In  contrast  to  a  window  location  closer  to  the  central 
maximum.) 

The  process  EXPOSURE  uses  to  determine  the  threshold  intensity  ratio 
can  be  called  “tuning*.  The  Micron  Eye  analog  of  tuning  would  consist  of 
three  steps: 

(1)  Determine  the  approximate  diameter  for  the  fiber  to  be  measured 
using  an  optical  shearing  eyepiece. 

(2)  Determine  the  Micron  Eye  location  parameters  L  and  x  by  the 
methods  in  Chapter  5. 

(3)  Vary  exposure  to  obtain  a  defined  relationship  between  the  three 
points  on  the  Micron  Eye  photograph. 

As  an  example,  return  to  the  case  where  L  =  30  mm  and  x^  =  5mm. 
(x^  is  the  inboard  x  location  of  the  Micron  Eye  array.)  Figure  32  shows 

the  “tuned*  relationship  between  three  points  on  the  Micron  Eye  array. 
This  relationship  can  be  expressed  as  a  ratio  of  a  :  b  :  c.  The  theta 
locations  at  the  inboard  and  outboard  edges  of  the  Micron  Eye  are  known, 
e.g.,  0i(wr  =  .165  radians  and  0^^  =  .300  radians. 
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Q*mr  6i  62  63  6*Ur 

Figure  32.  Three  Point  Spacing  Ratio 


This  means  that  8tnrrr  <0,  <  02  <  03  <  0^^,  . 

To  find  0t ,  02 ,  end  6Z .  run  EXPOSURE  for  the  diameter  being  tuned. 

EXPOSURE  outputs  the  theta  values  'or  the  optimum  exposure.  (See  Figure 
33) 


The  theta  values  are  related  to  x  locations  by: 

xt  =  L  tan(0,) 

Xj  =  L  tan(02) 

X3  =  L  tan(03) 

where  a  =  x.  -  x. 

I  wnwfr 

b  =  X^X, 

C  =  Xj  -  Xj 
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INTENSITY  RATIO  *10' 


VII.  DISCUSSION  OF  RESULTS 


The  simulation  was  conducted  for  three  diameters: 
d,  =  5.803  jim 

=  7.254  urn 

dj  =  8.705  4m  (where  d1  and  dj  are  ±20%  dp 

For  each  6V  three  levels  of  error  were  introduced:  \%,  2%  and  5%. 
Figure  34  shows  the  three  intensity  profile  curves  for  the  dr 


SIMULATION  INTENSITY  PROFILES 


Figure  34.  Simulation  Intensity  Profiles 


The  singled  Micron  Eye  window  wos  determined  as  previously 
described  at  L  =  30mn  and  x^  =  5mm.  The  exposure  was  tuned  with 

respect  to  the  7.254jim  perfect  data.  EXPOSURE  recommended  the 
Threshold  Intensity  Ratio  of  .0000370.  This  Threshold  Intensity  Ratio 
remained  constant  throughout  the  simulation,  for  all  diameters. 

Before  the  simulation,  all  diameters  were  tested  with  ng  error  using 
the  7.254jim  Threshold  Intensity  Ratio.  DIAFIND  recovered  df,  dj  and  dj 

exactly  (i.e.,  to  the  three  digits  accuracy  of  the  original  diameters). 

The  simulation  originally  began  by  collecting  thirty  data  points  for 
each  diameter/error  combination.  Twenty  additional  points  were 
collected  (50  total)  to  produce  meaningful  histograms. 

Figures  35  through  43  are  histograms  depicting  the  results  of  the 
simulation.  In  general,  the  expectation  was  to  see  a  decrease  in 
resolution  as  more  error  was  introduced.  It  was  also  anticipated  the 
resolution  would  be  best  for  the  7.254pim  case  (for  all  values  of  error) 
since  the  exposure  was  "tuned"  for  this  diameter.  Finally,  it  was  hoped  the 
method  would  be  more  accurate  than  existing  laser  diffraction 
measurement  methods  0?  0.5$). 
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Figure  35.  5.803jim  Fiber  (IX  error) 
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Figure  36.  7.254jim  Fiber  (IS  error) 
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Figure  37.  8.705jim  Fiber  (IX  error) 
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5.803  jim  fiber  (2%  error) 
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Figure  38.  5.803}im  Fiber  (2%  error) 
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7.254  jim  fiber  (2%  error) 
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Figure  39.  7.254jim  Fiber  (2%  error) 
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Figure  40.  8.705jim  Fiber  (2%  error) 
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Figure  41.  5.803jim  Fiber  (5%  error) 
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Figure  42.  7.254jim  Fiber  (5%  error) 
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8.6814 


8.705  jim  fiber  (51  error) 
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Figure  43.  8.7054m  Fiber  (5%  error) 


A.  ACCURACY  VERSUS  RESOLUTION 


The  histrogrems  show  the  accuracy  and  the  resolution  of  the  method. 
The  accuracy  is  associated  with  the  largest  spike,  and  related  to  the  size 
of  the  interval  called  the  resolution.  The  resolution  is  the  ability  of  the 
routine  to  “see"  a  difference  between  two  different  fibers.  For  example,  if 
the  resolution  is  ,005jim  the  method  does  not  discriminate  between 
8.705pmn  and  8.700ym.  Table  2  shows  accuracy  versus  resolution  for  the 
results. 


Table  2.  ACCURACY  VERSUS  RESOLUTION 
TERROR  DIAMETERjim  RESOLUTION  ACCURACY  (Res/d*ct)  % 


1 

5.803 

.004 

.069 

7.254 

.0042 

.058 

8.705 

.0060 

.069 

2 

5.803 

.0048 

.083 

7.254 

.0054 

.074 

8.705 

.0049 

.056 

5 

5.803 

.0105 

.181 

7.254 

.0070 

.096 

8.705 

.0136 

.156 
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The  date  show  a  decrease  in  accuracy  and  resolution  as  error  is 
increased.  It  is  evident  in  most  cases  that  the  7.254jim  results  are 
better  because  it  was  the  "tuned"  diameter.  The  only  exception  is  the 
accuracy  for  the  8.705jim  fiber  (2%  error)  is  better  then  the  7.254jim 
fiber.  Overall,  the  largest  error  is  .18  percent  which  is  less  than  one  half 
of  the  error  associated  with  the  manual  laser  diffraction  method. 

B.  ERROR 

There  are  two  contributions  to  error  in  this  simulation. 

The  first  can  result  from  programming/calculation  errors.  Many  trial 
runs  of  the  software  were  tak*n  to  minimize  the  likelihood  of  this  kind  of 
error. 

The  next  type  of  error  is  the  digitizing  error  encountered  in  an  actual 
experiment.  This  error  is  simulated  based  on  the  maximum  intensity  of 
the  profile  curve,  in  the  region  of  interest,  to  apply  random  error  equally 
to  all  points.  Hod  the  random  error  been  based  on  each  point,  points  with 
less  intensity  would  have  less  error,  and  points  with  higher  intensity  more 
error,  which  is  inconsistent  with  the  physical  digitizing  process.  Using 
the  maximum  intensity  avoided  such  a  condition  but  it  cannot  be 
ascertained  that  this  is  the  optimal  representation  of  the  physical  system. 


The  results  of  the  simulation  are  encouraging.  The  largest  error  is 
.10*  or  less  than  one  half  of  the  manual  methods.  Because  the  programs 
were  carefully  developed  and  tested  it  is  unlikely  they  contri'  'ted  to  the 
error.  Also,  much  effort  was  directed  towards  accurate  simulation  of  the 
physical  system  so  that  the  results  would  reflect  what  can  be  expected 
from  the  actual  experimental  measurements. 

The  simulation  demonstrates  an  increase  in  accuracy  two  to  ten 
times  better  than  that  currently  possible  by  making  manual  measurements 
with  laser  diffraction.  The  method  also  lends  itself  to  automation  which 
makes  it  attractive  for  quality  control  purposes  and  research  for 
materials  development. 


IX.  RECOMMENDATIONS 


This  study  has  shown  the  feosibilty  of  computer  aided  diameter 
measurements.  There  are  many  directions  future  work  can  take. 

A.  SOFTWARE. 

More  development  and  testing  of  the  software  could  result  in 
increased  accuracy.  It  would  also  be  valuable  to  implement  the  software 
on  a  smell  computer  (like  the  Macintosh)  to  allow  real  time  processing  of 
actual  data. 

B.  HARDWARE. 

There  remain  some  hardware  considerations  which  must  be  resolved. 
The  biggest  of  these,  perhaps,  is  the  accurate  positioning  of  the  Micron  Eye 
array  in  the  direction  perpendicular  to  the  laser  beam.  In  an  effort  to 
increase  the  resolution  and  accuracy  of  the  problem,  the  system  may 
benefit  from  two  Micron  Eyes. 
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APPENDIX  A.  FRAUNHOFER  DIFFRACTION  THEORY 


The  following  is  o  brief  discussion  of  Fraunhofer  diffraction  theory 
with  emphasis  on  aspects  of  the  theory  which  relate  to  this  research. 

A.  FRAUNHOFER  DIFFRACTION 

Fraunhofer  diffraction  (Figure  44)  results  when  light  approaches  and 
leaves  an  obstacle  or  aperture  in  the  form  of  plane  wavefronts.  iRef.  4:p. 
176]  The  light  source  and  the  plane  of  observation  in  effect  are  at 
infinity.  A  collimated  laser  beam  is  ideally  capable  of  Fraunhofer 
diffraction  because  its  beam  consists  of  parallel  rays  advancing  in  phase. 


Caatral  Maxima  Marfaraaea  Waites 


Figure  44.  Typical  Single  Slit  Diffraction  Pattern 


B.  THE  CLASSICAL  SINGLE  SLIT  EXPERIMENT 

The  simplest  demonstration  of  Fraunhofer  diffraction  is  the  single 
slit  experiment.  (See  Figure  45)  Parallel,  co'  imated  light  passes  through 
a  slit  of  width  a.  The  diffraction  pattern  is  visible  on  a  screen  located  a 
listance  L  from  the  slit.  An  observer  at  point  P,  moving  across  the 
screen,  sees  a  succession  of  maximum  and  minimum  intensity  points. 
These  extrema  are  the  result  of  constructive  and  destructive  interference 
of  the  light.  For  example,  a  minimum  occurs  when  the  angle  0  produces  a 
phase  difference  of  one  wavelength  between  the  rays  at  the  upper  and 
lower  edges  of  the  slit.  Thus,  minima  occur  whenever 

a  sin(6)  =  mA,  where  m=  1,2,  3 .  (A.!) 

These  minima  are  referred  to  as  interference  nodes.  (For  further 
discussion  of  this  subject,  see  Meyer-Arendt,  (Ref.  4].) 
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L 

Figure  43.  Single  Slit  Path  Difference  Relation 
to  the  interference  Node 


1.  Diffraction  Minima 

Examination  of  equation  (A.l)  shows  that  as  the  slit  becomes 
narrower  the  angle  0  becomes  larger.  As  this  theory  is  extended  to 
approximate  the  diffraction  phenomena  of  an  obstruction  (a  fiber),  one  can 
expect  short  distances  between  interference  nodes  for  larger  fibers  and 
greater  distances  between  nodes  for  fibers  of  smaller  diameters. 


Another  important  point  concerns  the  distance  between  minima.  It 
appears  the  interference  nodes  are  equidistant.  This  is  true  only  within 
the  limits  of  the  small  angle  approximation.  The  distances  between 
interference  nodes  actually  increases  as  one  moves  outward  from  the 
central  maximum  by  the  relation: 

=  Sin  l  m^/,d  1 

2.  Diffraction  Maxima 

The  diffraction  maxima  are  not  located  midway  between  minima. 
The  locations  of  the  maxima  can  be  derived  as  by  Meyer-Arendt  [Ref.  7:p. 
220].  These  occur  whenever  the  derivative  of  the  intensity  is  equal  to 
2ero: 


dl/dfi  =  dl/dfi  [  i#  (sin  B/B)2  1  =  0 

=  2I«  slnB/B  I  -sinB/B2  *  cosB/B  ] 


or  whenever: 


B  =  tanB 


(A.3) 


y=tanft  end  y=B.  (See  Figure  46)  The  maxima  are  displaced  slightly 
from  center  towards  the  central  maximum.  Much  further  away  from  the 
central  maximum,  the  maxima  are  nearly  halfway  between  minima. 


Figure  46.  Locations  of  the  Diffraction  Maxima 

C.  FRAUNHOFER  DIFFRACTION  AND  THE  FOURIER  TRANSFORM 

The  foregoing  discussion  centered  on  the  Fraunhofer  diffraction  due 
to  a  slit.  The  mathematics  of  the  slit  example  are  simple  and  provide 
insight  into  diffraction  physics.  A  more  elegant  derivation  of  Fraunhofer 
diffraction  shows  that  the  diffraction  pattern  of  an  object  is  the  Fourier 
transform  of  that  object.  (Ref.  9:p.  1741 


-b 


A 

Figure  47.  Fourier  Transform  for  the  Single  Slit 

Thus,  F(0)  describes  the  amplitude  of  the  diffraction  pattern,  and 
IF(0)12  represents  the  intensity. 

Modeling  the  transform  for  the  obstruction  is  more  complicated, 
consider  the  complement  of  the  slit  transform:  g(x)  =  1  -  f(x). 


79 


This  results  in  the  form  l  ••  -  2sin  (b 0)/0  ]  which  is  not 
transformable. 

An  alternative  is  to  use  a  substitute  function  h(x)  =  g(x)  -  f(x). 
Physically,  this  is  approximating  the  obstruction  as  two  parallel  slits: 


Figure  49.  Fourier  Transform  for  Two  Parallel  Slits 


which  has  the  solution: 

H(0)  =  (2sin  c0)/0  -  (2sinb0)/0 
where  c  would  be  large,  but  not  infinite. 

This  Is  still  a  crude  approximation  which  shall  be  Improved  upon  in 
the  next  section. 
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D.  REJECTION  OF  THE  SLIT  APPROXIMATION 

In  their  paper  titled  Tiber  Diameter  Measurement  by  Laser 
Diffraction"  [Ref.  3:p.  1378],  Perry,  Ineichen,  and  Eliasson  conclude  that 
the  diffraction  pattern  of  a  real  fiber  is  sufficiently  different  from  that 
of  a  slit  to  warrant  the  slit  approximation  being  treated  with  caution. 

Further,  they  recommend  a  solution  presented  by  Kerker  [Ref.5:p.  260] 
which  has  been  adopted  in  thi?  study.  Kerker’s  solution  assumes  the  fiber 
is  perfectly  reflecting  (i.e.,the  reflective  index  m=oo),  and  while  this  is 
not  completely  true.  Perry,  et  al.,  [Ref.  3:p.  1378]  indicate  some  degree  of 
absorption  is  not  likely  to  be  significant. 

E.  INTENSITY  EQUATION  FOR  A  REAL  FIBER 

Kerker  [Ref.  5:p.  260]  gives  the  scattered  intensity  relation  for  a  real 
fiber 

1/l#  =  (2/K,Ltf)  |  b,  *  2  I  b,  cos(ne)  F  (A.4) 

where  0  =  the  scattering  angle 

K#  =  2if  /X  (  X  =  laser  wavelength) 

ba  -  J,<«)  / 
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and 


a-  ifdj/A  (df  =  fiber  diameter) 

JB(or)  ar^  Bessel  functions  of  the  first  kind, 
to&Kae)  are  Hankel  functions  of  the  second  kind. 

The  real  fiber  equation  (A.4)  is  somewhat  obscure  in  Its  compact 
form.  It  can  be  shown  that: 

I/I  =(2/KLtf)  {(  r  *22  rcos(n0)  )2+(s +22  s  cos(n0)  J2]  (A.5) 

w  w  9  9  9  9 

where  rB  =  JB2(ar)  /  i  *  Yb2(qt)  1  (A.6) 

and  sB  =  Yt(ar)  /  [  J  *(cx)  ♦  ]  (A.7) 

Now,  the  intensity  ratios  for  any  0  location  can  be  calculated  for  any 
diameter  f1ber. 

1.  Sensitivity  of  the  Results  to  the  Number  of  Bessel  Terms 

The  calculation  of  Intensity  Ratios  requires  the  computation  of  n 
and  s#  Bessel  terms,  equations  (A.6)  and  (A.7).  Here  the  phrase 

"Bessel  terms"  indicates  the  algebraic  combinations  of  the  JB  and  Ya 
Bessel  functions. 
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o.  Two  competing  phenomena 

There  exist  two  competing  phenomena  which  govern  the  number 
of  Bessel  terms  to  be  used  in  the  calculations.  The  first  requires  a 
minimum  number  of  terms  for  accuracy.  The  second,  limits  the  number  of 
terms  so  the  Y  functions  do  not  cause  an  underflow  error  during 
computation. 

(1)  Minimum  Number  of  Terms.  As  in  any  series,  there  is  a 
minimum  number  of  terms  required  for  computational  accuracy.  Figure  50 
shows  the  effect  of  the  number  of  Bessel  terms  computed  for  seven 
curves,  all  with  a  diameter  of  eight  microns. 

Note  that  the  43,  50,  75,  and  86  term  curves  are  nearly  identical  and 
the  curves  with  fewer  than  43  terms  are  decreasing  towards  a  smooth 
line.  This  research,  indicates  that  50  terms  returns  four  digit  accuracy 
for  diameters  ranging  from  5  to  10  microns. 

(2)  Maximum  Number  of  Terms.  The  maximum  number  of  terms 

depends  on  the  value  of  oc.  As  the  diameter  decreases,  so  does  a.  This 
in  turn  produces  very  large  values  of  (in  the  rm  and 

denominator)  which  causes  an  underflow  error.  For  the  diameter  range 
considered  in  this  research,  underflow  was  not  a  problem.  In  one 
instance,  below  5  microns,  it  was  found  le  number  of  Bessel  terms  had  to 
be  reduced  to  43  to  prevent  underflow.  The  effects  of  the  minimum 
number  of  terms  was  not  investigated  at  this  smaller  diameter  range. 


INTENSITY  PROFILE 


Figure  50.  Effect  of  Number  of  Bessel  Terms 

2.  Effect  of  Different  Diameters 

Recoil  thot  the  diffraction  pattern  for  a  slit  showed  that  as  the 
diameter  decreased,  the  interference  nodes  moved  away  from  the  center  of 
the  pattern.  It  is  interesting  to  note  that  the  pattern  for  a  real  fiber 
behoves  the  some  way.  (See  Figure  5 1 ) 
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APPENDIX  B.  COMPUTER  PROGRAMS 

The  following  is  c  list  of  programs  discussed  in  this  Appendix: 

DATAMAKR 

DIAFIND 

EXPOSURE 

The  programs  are  written  in  Waterloo  Fortran  IV  (WATFIV)  and  run  on 
the  Naval  Postgraduate  School's  IBM  360  computer. 

A.  FORMAT  OF  DATA 

All  data  is  formatted  using  exponential  notation.  The  most 
frequently  manipulated  files  are  those  containing  theta  locations  and 
intensity  ratios.  The  format  for  these  files  is:  (IX,  E17.10,  IX,  E17.10). 
In  WATFIV  all  data  files  must  be  of  filetype  “WATFIV". 

To  compile  and  run  a  program  on  the  IBM  360  type  "  WATFIV 
PROGNAME  DATAFILE  *(XTYPE\  where  PROGNAME  is  the  filename  of  the 
program  and  DATAFILE  is  the  filename  of  the  data  file. 


B.  COMPUTER  PROGRAMS 

This  appendix  presents  the  various  computer  programs  and  outlines 
their  logic. 

1.  DATAMAKR 

DATAMAKR  generates  the  data  used  in  the  simulation.  The  equation 
for  computing  the  Intensity  Ratios  of  the  diffraction  pattern  for  any  value 
of  theta  is: 

l/l#  =  (2/K#Ltf)  |  b#  ♦  2  2  bB  cos(n0)  F  (A.4) 

where  0=  the  scattering  angle 

K#  =  2if  /A  ( A  =  laser  wavelength) 
b„  s  JB(or)  /  W&Hoc) 

and  oc-  udj/A  (df  =  fiber  diameter) 

J  (or)  are  Bessel  functions  of  the  first  kind, 

II 

H  <2>(ct)  are  Hankel  functions  of  the  second  kind. 

n 

All  computations  begin  by  calculating  the  required  number  of 
Bessel  terms.  Two  NONINSL  subroutines  are  called:  BESJ  for  the  J  (oc)  and 

n 
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BESY  for  the  Y JLce).  Note:  should  double  precision  be  desired,  the  NGNIMSL 


subroutine  BINT  will  return  double  precision  values  for  the  Jm  and  Yr 

The  Bessel  function  values  are  stored  in  an  vector  J(l)  and  Y(l), 
and  are  used  to  calculate  values  of  R(I)  and  SO).  Next,  the  matrix  of 
cosines  is  contructed.  The  size  of  this  matrix  is  M  by  K,  where  N  =  M-1 
and  N  is  the  number  of  Bessel  terms.  K  is  the  number  of  theta  locations. 
Values  of  theta  at  each  location  have  previously  been  stored  in  an  array 
called  TOO. 

To  visualize  the  program's  calculations,  the  matrices  symbolic  of 
these  calculations  are  sketched.  (See  Figure  52)  Note  that  r#  and  so  are 

outside  the  summation  term  in  equation  (A.1)  and  therefore  will  not  be 
multiplied  with  any  cosine  terms.  Since  array  subscripts  must  be 
designated  with  a  nonzero  value,  r#  and  s#  will  be  represented  by  r(l) 
and  s(  1 ). 

After  multiplying  the  two  matrices,  it  is  necessary  to  complete 
the  summation  by  summing  the  columns  in  the  matrix.  This  results  in  snN 
term  value  for  each  value  of  theta.  A  straightforward  calculation  then 
yields  the  intensity  ratios.  The  intensity  ratios  are  stored  in  an  array 
with  their  corresponding  theta  locations. 
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Figure  52.  DATAMAKR  Matrix  Algebra 


The  logic  of  DATAMAKR  is: 


Input - >  *  Bessel  terms  to  be  computed 

diameter  of  the  fiber 
screen  to  fiber  distance,  L 
array  of  theta  location  values 

Calculate:  Rn  and  Sn  using  BESJ  and  BESY 

Produce  the  Matrix  of  Cosines 
Produce  the  Matrix  of  Products 
Sum  the  columns  of  the  Product  Matrix 
Compute  and  output  the  Intensity  Ratios 

2.  DIAFIND 

This  program  finds  the  diameter  of  a  fiber  through  an  iterative 
process  of  residual  comparison.  The  program  accepts  the  data  which  has 
been  output  by  the  program  EXPOSURE.  The  user  is  prompted  for  K  and  df 

K  is  the  number  of  theta  values  the  program  should  expect  and  di  is  the 
initial  guess  of  the  fiber  diameter  (in  microns). 
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The  program  first  calculates  intensity  ratios  at  the  diameter  d., 
corresponding  to  the  input  theta  locations.  The  difference  between  the 
input  intensity  ratios  and  those  associated  with  df  is  called  the  residual. 
The  program  calculates  another  set  of  intensity  ratios  based  on  a  new 
diameter  d?  +  M.  where  Ad.  is  a  small  increment  of  diameter  (usually 

.05  microns  to  start).  The  intensity  ratios  calculated  using  dj  +  Ad,  are 

compared  with  the  input  intensity  ratios  to  give  a  second  residual.  The 
residuals  are  compared  and  program  logic  determines  whether  or  not  the 
Ad  increment  is  producing  convergence  to  the  actual  diameter.  The 
initial  guess  diameter  is  incremented  and  decremented  as  necessary  until 
a  desired  level  of  accuracy  is  achieved. 

A  key  to  understanding  the  convergence  process  is  the  residual  curve 
in  Figure  53.  It  is  important  that  the  initial  guess,  dj,  be  fairly  close  to 
the  actual  diameter. 

For  example,  the  residual  curve  in  Figure  53  is  for  an  actual  diameter 
of  7  microns.  If  di  is  greater  than  9  microns,  the  program  will  not 

converge  to  the  correct  diameter,  but  rather  to  a  diameter  just  above  10 
microns. 

It  is  hard  to  define  exactly  the  limits  of  di  for  any  given  actual 
diameter.  Figure  54  is  a  residual  cun/e  for  5  microns  and  shows  an  upper 
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limit  of  9  microns  for  d,.  A  thorough  study  to  define  limits  for  d4  has 
not  been  conducted,  although  convergence  has  always  been  attained  by 
guessing  d.  within  ±1  micron  for  diameters  ranging  from  5  to  9  microns. 


RESIDUAL  VS  DIAMETER 


Figure  54.  Residual  Versus  Oiemeter  for  5um  Fiber 


3.  EXPOSURE 

EXPOSURE  simulates  the  operation  of  the  Micron  Eye  by  taking  a 
photograph  of  intensity  ratio  data.  EXPOSURE  reads  in  the  perfect  data 
generated  by  DATAMAKR  and  digitizes  the  data  with  respect  to  a  threshold 
intensity  ratio. 

The  optimum  intensity  ratio  for  an  exposure  is  termed  the 
threshold  intensity  ratio.  The  threshold  intensity  ratio  is  located 
where  the  absolute  value  of  the  derivative  of  the  intensity  ratio  with 
respect  to  theta  is  a  maximum.  This  assures  that  the  threshold  intensity 
ratio  is  located  in  a  region  of  the  curve  which  is  closest  to  a  straight  line. 
This  reduces  the  effects  of  subsequent  interpolation  errors. 

The  threshold  intensity  ratio  is  calculated  by  the  subroutine 
DER1V.  Since  the  input  intensity  profile  curve  will  have  at  least  three 
maximum  derivative  points,  DERIV  calculates  the  average  of  the  three 
intensity  ratios.  This  average  is  then  passed  to  the  calling  program, 
EXPOSURE. 

EXPOSURE  searches  the  intensity  ratios  and  compares  them  with 
the  value  of  the  threshold  intensity  ratio.  This  process  corresponds  to  the 
Micron  Eye  addressing  a  pixel  and  comparing  its  voltage  to  the  threshold 
voltage.  If  an  intensity  ratio  is  greater  than  the  threshold  value,  the  ratio 
is  assigned  a  digital  value  of  1.  If  the  intensity  ratio  is  less  than  the 
threshold  value,  the  ratio  is  assigned  a  digital  value  of  0. 
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The  digitized  data  is  then  searched  to  find  the  theta  locations  at 
which  the  digital  intensities  change  from  0  -->  1,  or  from  1  — >  0.  The 
theta  locations  are  averaged  to  produce  a  theta  location  which  is  very 
close  to  the  threshold  intensity.  Because  this  averaging  process  is 
occuring  in  the  regions  of  steepest  slope,  error  is  assumed  to  be 
minimized  since  the  curve  can  be  approximated  as  a  straight  line  (refer  to 
Figure  30  in  Chapter  6). 

EXPOSURE  provides  the  user  with  the  option  of  introducing  error 
into  the  perfect  data.  The  program  prompts  the  user  to  input  the  desired 
error  and  calls  the  subroutine  RANDOM.  RANDOM  introduces  random  error 
into  the  intensity  ratios  and  returns,  the  imperfect  data  to  EXPOSURE. 
EXPOSURE  digitizes  the  imperfect  data  and  searches  for  the  occurences 
where  0— >1  and  1— >0.  imperfect  data  will  have  many  more  occurences 
of  the  digital  intensities  changing  from  1— >0  and  0— >1.  EXPOSURE  will 
average  the  theta  locations  if  they  are  in  the  same  location  (i.e.,  within 
+.005  radians).  This  method  has  been  tested  for  random  error  up  to  5 
percent. 
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CALCULATE  THE  R(N)  AND  S(N) 
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-*-x 

to 

• 

• 

Ol 

- 

• 

rH 

• 

W 

w 

- 

• 

x 

0$ 

EH 

• 

z 

- 

w 

W 

n 

W 

EH 

W 

w 

M 

W 

EH 

z 

Z 

o 

< 

W 

Z 

Eh 

z 

— 

O 

< 

w 

X 

n 

Q 

EH 

to 

5S 

W 

, 

Z 

S 

o? 

Z 

Eh 

i — ! 

M 

t-l 

z 

w 

P 

£ 

M 

X 

*"*x 

< 

X 

% 

•» 

P 

Eh 

a! 

P 

EH 

P 

w 

W 

p 

w 

W 

EH 

EH 

Z 

hH 

Eh 

P 

W 

< 

Z 

Jg 

t-H 

HH 

w 

X 

o 

o 

S 

M 

o 

Z 

2 

P 

w 

W 

<; 

x 

w 

EH 

Eh 

X 

W 

< 

W 

W 

Z 

W 

•  • 

Jg 

M 

X 

W 

w 

z 

Z 

P 

w 

EH 

D 

w 

w 

M 

w 

W 

O 

EH 

w 

P 

Z 

o 

2 

W 

z 

W 

EH 

<« 

P 

a 

w 

m 

Z 

n 

< 

w 

w 

H-t 

M 

P 

M 

n 

w 

W 

EH 

P 

w 

p 

X 

l-H 

z 

o 

o 

EH 

f*—** 

z 

P 

EH 

EH 

O 

m. 

HH 

< 

-  P 

< 

P 

O 

< 

Z 

EH 

2 

EH 

EH 

EH 

Z  M 

W 

Z 

P 

P 

O 

<— X 

S-Hh 

W 

n 

i — i 

b 

Z 

< 

tQ 

OrHM 

z 

W 

w 

z 

O 

EH 

WwZ 

o 

X 

W 

Pd  - — '  •— i 

EH 

1 

p 

ZWM 

w 

to 

Z 

z 

p  - 

- 

P-' 

P 

Z 

PH- 

- 

o 

Z 

a 

X 

Z“ 

rH 

z 

OP  ' 

p  ' 

1 

z  . 

t-H 

-  X 

X 

w 

X 

-"'cH 

r~M  ,p 

< 

r~rH 

too 

P 

wo? 

- 

to-  .< 

*  l-H 

o 

WO  ' 

O'  ' 

OO? 

M 

2H- 

rH— 

H  •>  .D 

P 

z 

SH- 

rH- 

.  1 

^-o 

o 

-  - 

O  rH 

x< 

X  .eh 

< 

a 

O  \ 

\ 

w 

rHMOS 

rH 

OOcn 

EH  - 

—  M 

M|| 

w 

EH  - 

O  II  rH 

-~WZ 

n 

rH  cn 

-  -  - 

ZtQ-wtOwP 

tO-^W< 

P 

ztOw 

tO-w 

. 

HH 

ii 

11  II 

v^EHZ 

PtwEH 

-~EH 

OP  II 

HZ  II 

p 

II  II 

XXX 

ftWWC  ' 

WCO  - 

II  M 

W< 

W< 

< 

P 

WEH 

zz 

EH  EH  EH 

PE H2 

HSQ 

EnSZP 

PEhS 

EhS 

II  EH< 

P-~- 

PZ 

OWWEH 

ZZZ 

Jh2 

MC&  < 

mWO< 

<Cd 

pm2 

MPi 

PEH 

IIPrH 

SprICSHW 

HHH 

<zozow 

ttiOHW 

MW 

<zo 

zo 

ZWW 

*C 

WO  11  POZW 

ZZZ 

USWSWZ 

ZWSjZ 

PP 

U2H 

SH 

ZPP 

POZ 

2QMWPPEH 

ZZZ 

oq 

LO 

tO 

O' 

rH  r- 

rH 

rH 

rH 

rH 

rH  rH  rH 

o  o  QUO  oo  oo  ooo  oo  o 


1G6 


WHILE(  ABS( TEST)  .  GT.  BETA)  DO 
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1-1 


I  F(  LOWER  .  EO.  0)  THEN 


Z 

W 

D3 

H 


Z 

w 

a 

EH 

rH  O 

o 

rH 

II  II 

II 

o  o 

cx 

o 

<  < 

w 

< 

j  a 

• 

w 

fa  fa 

EH 

fa 

ID 

z 

fa  CO 

D 

M 

O 

w 

o 

O 

cn 

QEH 

w 

So 

hi 

w 

WO 

n 

CS 

Z 

• 

z 

z 

w  . 

z 

w 

W 

W 

a 

w 

o 

a 

a 

EH 

a 

• 

EH 

EH 

EH 

w 

ID 

ID 

w 

ID 

ID 

O 

O 

o 

O 

1 

| 

( 

1  ID 

WiD 

W 

Q»H 

WiD 

WO 

mo 

rH 

5  1  "7 

mo 

H  l 

O  1 

O 

rfjMrH 

O  1 

OW 

OW 

o 

.wl 

.w 

.in 

•  rH 

♦ 

Wh 

rH 

o 

o 

^Ww 

.o 

.o 

.O 

• 

rHQarHO 

CX  .I-H 

CX  .rH 

CX  .rH 

a 

! — i 

W 

w 

W 

w 

.  II  II  II  II 

.11  II 

.11 II 

.11  II 

. 

+ 

CX 

CO 

CO 

a 

OJ 

w-—+- 

H 

QQQco 

QQJco 

QQtJcO 

9r 

EH 

•  rHrH 

J 

<<w 

«w 

<<w 

<.  Eh 

z 

w 

EhEhQO 

EHEHQO 

EhEhQO 

Ehh 

Q 

Q 

iJQZEh 

JWSEh 

JJSEh 

JD 

O 

ww 

3: 

WWWO 

WWWO 

WWWO 

wa 

rHO 

QQ 

W 

QQZO 

QQZO 

QQZO 

Q 

12 

Z 

"w* 

+  II 

W 

fa 

fa 

fa 

fa 

fa 

2 

fa 

m 

n 

i— i 

n 

l-H 

i-H  Eh 

_ _ * 

M 

W 

W 

W 

w 

Z 

U1 

CD 

W 

w 

Q 

ii  a 

— ✓ 

Q 

a 

Q 

Q 

tJ 

Z 

o 

hi 

Z 

w 

W 

W 

w 

W 

MO 

HH 

W 

in 

co  to 

O  O  DO  O  UO  O 


OO 


O 


o 


lOCO 

coco 


OO 


o 


OO 
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rHrHO 
.11  II 

a 

WH 

.  j 
w 

j  p 

w  2 

gu 

z 

u 

z  t 


Q  — 
Z  fa 


2  H 

z 

Q  m 
Z  Pi 
fa  fa 


*  * 

*  * 

<  * 

*  * 

*  * 

*  * 

*  * 

*  * 

*  * 

*  * 

*  •* 

*  * 

*  * 

*  * 

*  * 

*  * 

*  * 

*  * 

*  * 

*  * 

*  * 

+  * 

*  * 

*  * 

*  * 

*  *  fa 

*  -it  J 

*  *  P 

*  *  a 

*  *  o 

*  *  s 

*  * 

*  *  H 

*  *  P 

*  *  fa 

*  *  H 

*  *  P 

*  *  O 

*  * 

*  * 

*  * 

*  * 

*  * 

*  * 

*  * 

*  *  i 

*  * 

*  * 

*  * 

*  * 

*  * 

*  *  i 

*  it 

*  -it  ( 

•e  *  < 

*  *  I 

*  * 

*  -it 

*  -it 

*  -it  i 

*  * 

*  *  i 

*  * 

*  -it 

*  -it 

*  * 

*  * 

*  -it 

oooooooo 


<  • 
*  H- 


H  iJ  Z 

fa  <  O 

£  M  M 

<  H  H 

hH  HH  «C 

Q  Z  fa 


fall  W 
fa  fa 
H 

<n  u 

Hw  ^sn 


fa 

'■Sd  <  U 


fcd  faH  hH  fa 
fa-  — -  fa¬ 
ct!  H  '  fa  '  Q 
O  -  -  -  -  ■ 

fa  '  ' 


Id  hS  P  EG  , 
{C  -*-H  H  H  fa 
H  fa  0<  < 

-  Id-  <►-*-  H 
-  Q  '  Q  '  O 

-  -  -  -  fa  -  -  H- 

-  EC  - 


o  oo 

<  HU 


Id 

w 

Id 

fa 

fa 

J 

. 

E 

tJ 

U 

fa 

p 

fa 

H 

< 

H 

Id 

Q 

p 

U 

H 

P 

o 

II 

N 

fa 

H 

£ 

Z 

O 

g 

rH 

U 

O 

4, 

tD 

< 

w 

H 

< 

HH 

o 

Q 

H 

H 

Q 

o 

P 

fa 

fa 

CO 

J 

P 

fa 

o 

EH 

II 

O 

fa 

fa 

- 

fa 

a 

•J 

n 

w 

Cd 

o 

fa 

< 

p 

II 

N 

% 

hH 

o 

fa 

o 

cx  faH 

“  .S 

fa  uu 

n  OH 


fa  Z< 
H  PH 


-  -  O —  t-h-  -  C'l-  -  Hro- - O' 

faw  fa-w-  fa—  -  fa—  O' 

„  ,  h  .  Eh  >  H  '  '  H  '  '  "O' 
HHH<HH<HH<HHH<HHH 

zzzszzszzszzzszzzp 

hH  hH  t-l  OS  M  HH  05  fa  *-<  *-<  fa  M  *-*  H 

dosoiofafaofafaofafafaofafafao 
pH  Cu  Hi  fa  fa  fa  fa  fa  fa  fa  fa  fa  fa  fa  fa  fa  fa  O 


oc 

-  O'  OH  fa 

O'  -  -  P  fa 

'O'  '  '  z  ^ 

H  HH  *h  E 

ZO  ZZ  H  H 


HhEH  I — I  < — I  Z 

OSO  fafa  o 

PhO  fafa  o 


O' 

t"-  O  O' 

cn  ^  O' 

uoa o  oo  oo 
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WRITE  STATEMENT  WRITES  THE  RESULTS  TO  THE  FILE 


XX  ' 

f-irH 

Vx  CO-  O 


Q  O' 
rHO 


sX  co-  o 

'  H  Bi 

5  '  W 
HO  'll 
■rt  —  CO 


Eh  t"  .  U1U  —* 
Z  iHr^  it  SO  O 
P  WrH  Sz  r-t 

O  w  mw^ 
a  '  zh<2  t" 

-  *  o  rH 

*  -  *Hj<  W 

II  II  HWH  ' 

*-i  <cow  X 
_  <05  o5coa  h 

W  HW  ■WWH  ' 
W  WHXH«-  - 

Q  COWrHt-H 

£  -  II 

•*-<'-  rH 

<  l-l\  -  ^ 

»H  "OXX  '  '  rH 
Q  X  r-lX\  + 
•^■rHD  '  «rH\ 

<  +  'OH\V)  W 
pn\H  .  '\n  W 
EHw  Ot^O  '  O 
OWO<rHrHO  - 
<tdr+-  W  .rH-  t0“ 

<5  •  '  t''  •  z 

r"  t^X  '•HI''  <  ' 
HH  mHi-H-  WrH  II  HO 
^_^Sd  W  QrH 

05H  'll  '  «  • 
+  '\  -  'HWr- 
HH-  Xa5  -  W  rH 
<W  '4x3  II  a  w 
EhP5  II  OEh  II  EH-  ' 
W  rHW05  -X 
CQ  -H  .SUQtw^'rH 
r~C0C"<H  O  •  ' 

>XCdrHHH[d<  rH- 

EH— HWDSHCOrH 
COE-t-  <JMW  II 
W  '  -COhhWP  ' 

Eh  C0QQD-  — 
- — <H  W  <  rH 

r'rH  'I|DQiJ>0  + 
t^'-r  Ofd<  Eh*-* 
i/>HXJ  HZ# 

'  <  J<HH-  -  05 

oo  'wP<JW  '  ' 
-^XEhQhhP-  X'tf- 

W  <HH£hO  rH  .  ' 

£H  '  -rH\ 

n£05WZ<XXrH\ 
05  OC5hOH\W  ' 
IE  4w-  -  - 
US  USUSUS'X'X'X'X 


US 

o 

HHtH 

H  CO 

< 

OCXQ 
H  u 
>*Z 
Q£h2 

W-O 
Ed  CO  (OEH 
JZPZW 
(OEH  Cd05 
S  tdEH 
hhWOSZW 
ZH<nnp 

0<  J 
ZJCOO< 
DZS> 
CdOOEH 

05J*-i  W 

•<<eh  a 

Q  OO  -EH 
M>H  ZCO 

cq  cncoaboco 

td  WWOWhhi-h 

hh  aw<a 

EH  PQStoSo 
b  CdZ  Cd  m 

o  a<-atHa 
a  j  zhhS 
a  <ao  w 

P  OCOhh  COiJ 
CO  Mto  .Z< 
cqwcowp 

£  US  OZHQ 
<  rH  [dOZH 
05  OSrHMCO 

O  IE  .(OEh  W 
O  OS*  OW05 
05  dcnuza 

Ph  W£dJbEH< 

«mow 

Z  WW 
QQhhJO> 
WZC0W  HH 

C0<'— cocoo  . 

a  cow  q 

acowbOHH 

cowcqdhco 
wwz  <  w 

Z«H[d>QOi 
hh  Ena  W 
HCOPHWCOW 

pjo  aaz 

OJ05WH  hh 
05<  O  ZEh 
CQO><  QWP 

a  oscozao 

C0Q<  Whh£ho5 

mojpw 
wen  cap  ww 
wwM<oo5a 

H05J>H<Eh 


J 
< 
cop 
CO  EH 

wo 
o5  a< 

W  .  O  ' 
H  W  <Q 
W  ZQ 
SCOhh  05 
<WEH05O 
hhhDOW 
DHOW 
rH  oi  COCO 
DCO  WW 

<zoaa 

pWZrU 

HEHM« 

OZJ» 

<HHp 

<>•>• 
W05UHH 
XW  l-HHH 
WHHOC0C0 
Z  WHZZ 
hhWS  WW 
HO<U5HH 
p  hOZZ 

OZQ<hh 

050  as 

CQhhCO  oo 

ancoozz 

C0<WWHH 

HPCOZZ 

oaoconHH 

Hw  «< 
SWWHH 

Doa  zz 

WOHWOO 
us  poo 

C0O05D 

<ZO<>^>H 

WHH 

H  '  «05 
05<C0D0505 
WHW<« 
HOhhP 

WHHEHDrHrH 

SQhh 

<  cocoxx 

hOZW 
Q<W05OO 
DH  OO 
WWZWrHrH 

a  HHP 

H<  H« 
coco  cococo 

HH  HHH 

C0H 

wo 

o< 

Qa  05HH 
05ZZ 


ooooooooooooooooooooooo  oooooooo 
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w 


/pass: 

TIONS 

csa 

OO 

oa 

> 

az 

OH 

HD 

•H 

O 

<a 

< 

O 

a 

><a 

H 

rH 

Da 

a 

— 

oa 

aa 

** 

s 

aw 

>a 

< 

W 

<w 

an 

►H 

ua 

a 

O 

w 

Q 

a 

w  _ 

030  Zh 

z  ta 
Oh  sa 
hx  aw 

OJCQ 
IOQ  On 

hz  za 
zo  n 
uo  a 
za  a3 


a 


op 

os 

zz 

o 

wa 

aa 

<H 

D 

HW 

aw 

aap<a 

Han  h 
«a 
aaans 

OODW 

UKZ 

aaaHn 

fOCOUwW 

ss  oo 

Zzo3><<! 

aaoeiw 

KKHan 
HH  <P 

win  incn 

t— «•— I  nn 


as  HJ 


<3 
>OH 
X  OrHW 

o  -a 

-  mHOO 

- - OH 

OHrHZ 

-  ow— n 
n  hhO  - 
a  '-QW- 

win  o 
<w  -  -o 
HO 

WH--00— 

zzoooa 

OhOhhh 
O  a—  <5 
SOH 
4HOJOOZ 
HU  ZHn 
o<  -a— 
ZH^OO  ' 
azo  a-* 

hO  '  O 
»  -O 

OJxOO- 
<3  riOH 
•st;  —hu 
q  *c 

aaaoa  -h 
««owoz 

S^  HhOh 

Z  ao"  — » - 

-  .*h  • — O- — • 

HEQZO  o 
<HOO  -O 

OSt-tZ*~,rH'~r-< 

apnw— o— 
o  z>hos 
UJJIjJBJhD 
H«sa— in 

ZWWhhhhHW 

HQ!0JQ 

<x<x<x 


w 


a 

< 

a* 

03 

EJ 

< 

a* 

\H 

O 

nX 

a— 

*  \ 

oo 

OQ 

OO 

cscs 

ii 

II 

H 

HW 

OZ 

ZO 

au 


a 

x 

H 

a 

o 

a 

H 

z 

a 

s 

o 

< 


a 


a 

a 

a 

Q 

a 

a 

H 


a 

x 

H 


< 

O 

03 


W 

a 

d 

a 

< 

> 

z 

o 

n 

B 


a 

a 

a 

w 

w 

a 

03 

O 

a 

a 


a  •  d 
3  aw  a 
a  <z  a 
o  a 
\  an 
oh  a 

—  UK 
Q  az  H 
DD 

♦  aa  a 
<  h 
>a  < 
a  a 
ww  a 
aw  z 
na  a 
a  o 
w 


x 


a 


a 

a 


Q 

Q 

CQ 

Z 

x' 


a 

a 

i— i 


ii  ooii 

n  WOO 

ho  a  —o 

O  I  03  II  n 
On  ||  —II 

n  a—o 
ii  &>  anaa 
o  a  <— aa 

QZn  QHjnn 


UOOUUO 


OO  UUUU  uu  uuou  uuoouuuo 


u 


u 
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« 

-it 

* 

W  Q  • 

* 

to  ZW 

* 

O  <05 

* 

z  .  w 

•It 

H— •  <H 

-It 

JC5W 

* 

o\ws 

* 

HHH 

*  OSmWCO 

*  CQOvD 

*  O—  in 

•It  HfSS  • 

*  S«CO 

*  mH05C/1 

*  otuemm 

*  05 

*  -<WK 

*  <11050 

■#  H<WH 

*  WHHZ 

*  swus 


*  Ehm  - 
■n  05  QW 

*  OWJO 
•*  twOSj-jZ 

*  U2< 

*  was  h 

*  wjom 

*  D  St-t 

*  J  -HQ 
■It  <H 

■It  >Z>ttC 
■n  W05H 

*  asw< 

*  wm>i4 

*  B5C5W 

•It  Mtd  O 

*  DOtWH 
•It  OXm 


•It 

05  WW 

X<S 

■It 

Z  05 

•It 

QOW 

•It 

s  w 

•H 

W<Z05 

O 

•It 

X  M 

* 

Z0<0 

•O 

It 

* 

HHHW 

T — 1  . 

•* 

>T><|I 

* 

ocn 

II  * 

•it 

WCQ 

W  -It 

W<wW 

hhW 

•It 

W  —  D 
CQ  II  *-•  Z 

w  M 

Jr-v*  h 

iJmp£  Z 
<— «  O 
OJHt-1  O 


•It  <U3  H 

•k  Ct'.  *T1  s 

•it  whww 
•It  Z  05m 

■It  wzw 

■It  OhHJ 


oiiz 

O  M 
<N-~H 
h3 
0—0 
QHO 


N 

* 


*  ^ 

■It  J 


Mw 

■it 

w 

■» 

* 

— s 

* 

z 

* 

HSO 

* 

i— t 

-*•* 

* 

— OZ 

* 

cn 

rH 

* 

z 

><zw 

* 

o 

1 

* 

— -* 

WQ 

* 

u 

t-H 

* 

cn 

Q 

■it 

Cn 

•— 

* 

+  \ 

* 

* 

•it 

o 

■ - ■ 

* 

Q 

■it 

X 

cn 

* 

* 

■it 

* 

§ 

•it 

M 

o 

* 

•it 

05 

o 

* 

<*"S 

. — v  w 

•it 

H 

* 

z 

WCNt« 

•n 

< 

II 

* 

— It 

* 

g 

* 

o’ 

>->■*  It 

* 

sx 

* 

* 

W 

* 

w 

* 

3 

■K 

rH  '7>h0 

It 

W 

II  II 

X 

II 

* 

H 

hZ 

rH 

W 

* 

< 

1 

»H 

Q 

* 

J 

OCn 

KH 

o— II  II 

Z 

•It 

D 

ocn 

— 

os 

t-4 

•It 

O 

tfcn 

u 

nO- — • 

H 

* 

J 

Zmw 

Z 

* 

< 

OO 

OW— 

O 

* 

O 

QQ 

QQ05W 

U 

•it 

j  ^  * 

OOUOUOUOOOOOOOOOOOOO 


O  OOOOOO  o  o 


•It 
•It 

1  * 

I  -It 
;  -it 

t  -It 
I  * 

;  -it 

>  -it 

i  * 

* 

* 

* 

>  -it 
)  * 
i  -it 

ouoo 
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GENERATE  THE  MATRIX  OF 


Er 2 

-Z  ' 

CS  M 
tl  CT>  I 

H(J\  M 

O  ■ — * 

oo  o 

rnQ  05 


■It 

CO 

■k 

Cd 

•It 

O 

■k 

V— 1 

•k 

05 

•k 

H 

z 

Z 

■k 

<J 

■k 

2 

% 

•k 

r-\ 

H 

* 

H 

1 

•k 

O 

M 

■k 

D 

•k 

Q 

* 

•k 

O 

O 

o 

■k 

05 

■k 

0i 

■k 

* 

* 

•k 

Cd 

•k 

EC 

* 

H 

M 

►H 

■k 

•k 

Z 

& 

CO 

•k 

n 

•k 

•k 

CO 

*  Z 

*  e 

*  P 

*  a 

+  o 

*  o 

* 
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